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ABSTRACT 


The  evaluation  of  the  beam  steering  capabilities  of  the  37-element  short-period 
seismograph  array  at  TFSO  indicates  that  the  array  provides  a  signal-to-noise 
ratio  improvement  over  a  single  detector  of  approximately  12  dB  at  frequencies 
near  1.0  cps .  The  phased  summed  output  of  the  array  suppresses  the  noise  by 
a  factor  ofV^  (15.6  dB)  for  frequencies  greater  than  about  1.5  cps.  At 
1.0  cps,  the  average  noise  reduction  relative  to  an  individual  output  is 
13.7  dB.  The  signal  amplitude  loss  in  a  30-second  P-wave  coda  is  between 
1.0  and  2.0  dB  in  the  frequency  range  of  greatest  signal  power  (0. 5-1.0  cps). 
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PRELIMINARY  EVALUATION  OF  BEAM-STEERING 
CAPABILITIES  OF  THE  TFSO  37-ELEMENT  ARRAY 


1.  INTRODUCTION 


In  November  1967,  construction  and  instrumentation  were  completed  on  a  medium 
aperture  37-element  short-period  seismograph  array  at  the  Tonto  Forest  Seismo- 
logical  Observatory  (TFSO).  The  design  of  the  array  was  directed  toward 
providing  the  capability  for  velocity  filtering  of  teleseismic  signals  and, 
more  specifically,  to  increasing  the  signal-to-noise  ratio  (S/N)  of  processed 
outputs  in  the  frequency  range  of  0.5  to  2.0  cps .  The  results  of  noise  studies 
conducted  on  data  recorded  by  the  cross-linear  and  31-element  arrays,  used  in 
conjunction  with  theoretical  array  response  studies,  indicated  that  significant 
S/N  improvement  could  be  provided  by  velocity  filtering  of  a  hexagonal  array  in 
which  the  distance  between  sensors  is  5  kilometers  and  the  overall  diameter  of 
the  hexagon  is  30  kilometers.  The  array  was  designed  on  this  basis;  however, 
due  to  considerable  topographic  relief  in  the  area  and  problems  of  legal  access, 
the  actual  locations  of  the  elements,  shown  in  figure  1,  were  slightly  perturbed 
from  the  ideal  theoretical  locations. 

This  report  deals  with  the  evaluation  of  the  beam  steering  capabilities 
exhibited  by  the  37-element  short-period  array  during  the  first  months  of 
its  operation. 


2.  DATA  SELECTION  AND  PREPARATION 


2.1  METHOD  OF  EVALUATION 

The  evaluation  consisted  basically  of  the  comparison  of  the  S/N  of  the  phased 
summed  outputs  of  the  37-element  array  and  of  two  subarrays  to  the  S/N  of  the 
average  of  the  outputs  of  the  individual  seismographs  and  to  the  unphased 
summation  and  the  unphased  filtered  summation  of  the  cross- linear  array. 

Figure  2  shows  the  locations  of  the  seismographs  which  were  contributing 
elements  to  the  unphased  summation  and  to  the  unphased  filtered  summation 
of  the  cross- linear  array.  One  of  the  subarrays  chosen  for  evaluation  con¬ 
sisted  of  the  center  seismograph  and  the  seismographs  of  the  inner  ring;  the 
other  subarray  was  composed  of  the  center  element  and  the  elements  of  the 
inner  and  middle  ring.  Twenty  earthquakes  occurring  at  teleseismic  distances 
from  TFSO  were  used. 

The  S/N  for  each  of  the  various  outputs  was  computed  by  two  methods.  The  first 
method  utilized  the  ratio  of  the  power  density  spectrum  of  the  signal  to  that 
of  the  noise;  the  second  used  the  ratio  of  an  amplitude  measurement  from  a 
dominant  half-cycle  of  the  signal  and  the  rms  of  the  noise.  Noise  reduction 
and  signal  loss  were  computed  in  addition  to  the  S/N. 
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Figure  2.  Cross-linear  array j  Tonto  Forest 
Seismological  Observatory 
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2.2  DATA  SELECTION 


The  data  samples  were  selected  from  16-millimeter  film  seismograms  of  the 
individual  elements  of  the  37-element  array.  The  daily  TFSO  station  message 
to  the  Coast  and  Geodetic  Survey  (C$GS)  and  the  daily  LASA  bulletin  were  used 
as  guides  in  the  search  for  suitable  data  samples.  Each  sample  selected 
satisfied  the  following  requirements: 


a.  At  least  5  elements  from  the  inner  ring,  10  elements  from  the  middle 
ring,  and  15  elements  from  the  outer  ring  were  operational  and  not  excessively 
noisy  at  the  time  the  signal  was  recorded.  Also,  it  was  required  that  each  of 
the  8  elements  of  the  cross-linear  array  was  operating  properly. 

b.  Each  data  sample  was  3  minutes  and  20  seconds  long,  consisting  of 
2  minutes  and  50  seconds  of  noise  followed  by  30  seconds  of  signal  coda. 

c.  Each  signal  originated  within  the  epicentral  distance  range  of 
30  to  90  degrees  from  TFSO. 


Many  samples  which  satisfied  these  requirements  were  later  rejected  when 
digital-to-analog  plots  of  the  Astrodata  digital  magnetic  tapes  indicated 
that  an  acceptable  number  of  channels  had  not  been  suitably  recorded.  Excessive 
spiking  and  parity  errors  on  the  digital  recordings  caused  many  samples  to  be 
rejected.  The  spiking  was  due  primarily  to  lightning  and  to  the  fact  that 
both  data  and  power  were  transmitted  over  the  same  spiral-4  cable.  The  parity 
errors  resulted  from  a  defective  tape  transport.  Ultimately,  20  data  samples 
satisfying  the  selection  criteria  were  obtained  and  processed.  Table  1  gives 
the  epicentral  data  for  each  signal  and  lists  the  unusable  channels  for  each 
sample.  The  epicentral  data  were  obtained  from  the  "Preliminary  Determination 
of  Epicenter"  cards  published  by  the  C$GS.  Figure  3  shows  the  distance  and 
azimuth  of  each  epicenter  relative  to  TFSO. 

The  calibration  constants  for  each  of  the  individual  outputs  and  for  the 
unphased  filtered  and  unfiltered  summations  were  determined  from  1.0  cps 
sinusoidal  calibrations.  Prior  to  24  May  1968,  the  individual  elements  were 
calibrated  both  separately  ard  simultaneously  in  groups  of  about  10  elements. 
After  24  May,  the  simultaneous  calibrations  were  discontinued  because  of  some 
doubt  as  to  their  validity.  The  calibration  constants  for  data  samples 
1  through  11  were  computed  from  the  simultaneous  calibrations,  while  the  cali¬ 
bration  constants  for  the  remaining  data  samples  were  necessarily  computed 
from  individual  calibrations.  As  a  check  on  the  two  types  of  calibrations, 

3  data  samples  were  processed  with  calibration  constants  computed  from  both 
individual  and  simultaneous  calibrations.  The  output  data  showed  that  the 
two  methods  of  calibration  gave  essentially  the  same  results. 

Digital-to-analog  plots  of  the  calibrations  were  obtained  to  determine  whether 
any  of  the  Astrodata  channels  were  overloaded.  In  some  cases,  overloading 
caused  the  calibrations  for  some  elements  to  be  recorded  improperly.  The 
problem  was  essentially  eliminated  by  the  installation  of  an  automatic  overload 
reset  circuit  in  the  Astrodata  system  on  16  April  1968. 
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Figure  3.  Location  of  signal  epicenters  relative  to  TFSO 
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Because  of  parity  errors,  the  calibration  constants  for  data  samples  1  through 
8,  excluding  sample  6,  were  computed  from  calibrations  recorded  on  a  day 
different  from  the  day  on  which  the  data  sample  was  recorded.  In  almost  all 
of  the  samples  there  were  a  few  elements  for  which  calibration  constants  were 
estimated.  Estimates  were  based  either  on  previous  calibrations  or  on  the 
basis  of  calibration  levels  of  other  systems. 

Calibrations  were  demultiplexed  from  the  Astrodata  field  tape  and  analyzed  by 
computer  to  obtain  the  calibration  constant. 

3.  DATA  PROCESSING 


After  the  calibration  constants  were  computed,  each  data  sample  was  demultiplexed 
from  the  Astrodata  digital  field  tape,  and  an  X-Y  plot  of  each  output  was  gene¬ 
rated.  Figures  4  through  8  show  a  portion  of  the  plots  for  each  individual 
element  and  for  the  unphased  filtered  summation  output  of  the  cross- linear  array 
from  data  sample  6  (see  table  1).  Each  trace  has  been  normalized  by  the  plot 
routine  on  the  basis  of  the  maximum  amplitude  for  that  particular  element; 
consequently,  no  relationship  between  relative  trace  amplitudes  can  be 
obtained  from  the  plots.  The  plots  were  used  to  determine  the  relative  delay 
times  necessary  to  form  the  phased  sums.  The  element  with  the  earliest  signal 
arrival  was  used  as  the  reference  element. 

Formation  of  the  phased  sums,  determination  of  signal  amplitude  and  calculation 
of  noise  rms,  calculation  of  the  average  output  of  the  individuals,  and  com¬ 
parison  of  the  average  output  with  the  phased  sums  were  accomplished  by  a 
computer  program  (TFOSAN) .  Because  three  phased  sums  were  desired  for  each 
sample,  that  is,  the  7-element  subarray  (PS1) ,  the  19-element  subarray  (PS2)  , 
and  the  37-element  subarray  (PS3) ;  the  sample  was  processed  three  times  by 
program  TFOSAN.  On  each  pass,  only  those  elements  contributing  to  the  parti¬ 
cular  phased  sum  under  consideration  were  used.  Besides  the  subset  tape 
containing  the  data  sample,  program  TFOSAN  required  as  input  the  delay  times, 
calibration  constants,  filter  coefficients,  label  of  any  element  to  be  deleted 
from  consideration,  time  window  length,  and  a  reference  start  time  for  t  e 
time  window.  The  filter  coefficients  used  were  for  the  standard  SDL  bandpass 
filter.  Figure  9  shows  the  amplitude  response  of  the  SDL  filter  and  the  analog 
filterused  in  the  cross-linear  filtered  summation  (CLSF) .  The  SDL  filter  was 
applied  both  forward  and  backward  in  time  in  order  to  achieve  zero  phase-shift, 
thus,  the  effective  bandpass  applied  to  the  data  is  the  square  of  the  SDL  filter 
amplitude  response.  Each  element  was  filtered  prior  to  the  ana-ysis  o  t 
data  sample.  The  phased  sum  was  formed  from  the  filtered  elements.  Figure  10 
shows  a  portion  of  the  plot  of  the  filtered  and  aligned  traces  for  Z1  through 
Z10  for  data  sample  6.  The  relative  amplitudes  between  traces  are  correct. 
Comparison  of  the  filtered  data  with  the  raw  traces  shows  that  the  signal  is 
enhanced  considerably  by  the  bandpass  filter. 

Program  TFOSAN  allowed  any  element  in  the  data  sample  to  be  deleted  from  the 
analysis  procedures  and  the  phased  sum.  This  option  was  used  in  forming 

PS1  and  PS2 . 
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Figure  4.  Calcomp  plot  of  Astrodata  recording  of  a  portion  of  data  sample 
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Figure  6.  Cal comp  plot  of  Astrodata  recording  of  a  portion  of  data  sample 


Figure  7.  Calcomp  plot  of  Astrodata  recording  of  a  portion  of  data  sample 


Table  1.  Epicentral  data  and  unusable  elements  for  each  data 


i/i 

QJ 

rH 

g4 

aJ 

</) 


m 

"3 

* 

CO 

im 

m 

*3- 

1^*. 

CM 

ro 

CM 

CM 

CM 

o 

ro 

ro 

ro 

o 

o 

ro 

no 

• 

C 

IM 

O 

M 

CM 

ro 

ro 

PM 

PM 

CM 

CM 

IM 

M 

no 

o 

Of 

CM 

M 

SJ 

M 

rM 

M 

CM 

CM 

E 

m 

* 

IM 

•> 

•* 

•» 

« 

0> 

» 

M 

IM 

QJ 

CM 

ro 

CM 

•* 

m 

m 

CM 

in 

• 

ro 

IM 

ro 

CM 

M 

CM 

CM 

■ 

o 

CM 

CM 

oo 

oo 

CM 

CM 

» 

QJ 

IM 

IM 

CM 

M 

IM 

p — - 

s 

ro 

M 

IM 

p— 

— 

rM 

M 

O 

CM 

r — • 

si 

M 

Ml 

IM 

CM 

r— * 

(U 

•* 

* 

IM 

#* 

0> 

m 

•* 

m 

» 

PM 

IM 

im  rM 

O 

CM 

m 

im  im 

•> 

o 

•» 

03 

ro 

m 

IM 

M 

IM 

ro 

O 

_Q 

r—  CO 

CM 

*  CO 

— 

, — 

ro 

ro 

^3" 

p — - 

CM 

—  ro 

t— 

r— 

ro 

no 

CM 

» 

•* 

ro 

IM 

05 

IM  M 

im 

03  IM 

M 

IM  M 

M 

t - 

CM 

M 

M 

—  IM 

r— 

IM 

M 

M 

M 

03 

in 

IM 

f— 

l/t 

IM 

M 

NJ 

M 

M 

M 

IM 

IM 

IM 

3 

m  m 

m 

• 

* 

m 

•> 

#• 

» 

* 

M 

m 

« 

m 

* 

C 

CM  CO 

CM 

*  ro 

p— 

r— 

r— 

^3- 

* 

» 

co 

«  VO 

m 

IM 

in 

IM 

03 

in 

IM 

» 

» 

ro 

0> 

ZD 

i—  ro 

r—  CM 

— 

p— ^ 

CM 

CO 

CM 

CO 

CO 

p—  ro 

lo  ro 

co  ro 

ro 

CM 

CM 

CO 

i^— 

CM 

CO 

M  M 

IM 

IM  IM 

IM 

M 

M  M 

I'M 

M 

M 

M 

M  M 

fM 

IM  IM 

IM  IM 

IM  IM 

IM 

IM 

IM 

IM 

IM 

IM 

IM 

o 

oo 

LL. 

o 

o 

0 

0 

D 

o 

o 

0 

o 

0 

o 

0 

D 

o 

o 

o 

o 

o 

D 

o 

mi 

03 

CM 

, — 

CO 

CM 

*3- 

^3- 

03 

p — 

oo 

oo 

p— 

p— 

X) 

00 

O 

*3- 

r— 

CM 

<  E 

, — 

p— 

03 

ro 

CM 

ro 

CM 

ro 

in 

ro 

CO 

CM 

^3- 

CM 

ro 

*— 

r— 1 

CM 

■ 

»“ 

ro 

ro 

•“ 

00 

ro 

»“ 

CM 

CM 

r“ 

*“ 

r— 

ro 

ro 

ro 

ro 

ro 

4- 

o 

oo 

•  u_ 

o 

o 

0 

o 

s 

o 

0 

0 

0 

o 

o 

3 

0 

0 

O 

D 

0 

o 

3 

0 

■M  t— 

CO 

O 

o 

oo 

LD 

ro 

r— 

^3- 

CO 

CO 

oo 

O 

oo 

CM 

03 

ro 

PM 

oo 

l/) 

« 

« 

• 

• 

• 

• 

■ 

• 

• 

• 

• 

■ 

• 

• 

• 

• 

• 

■ 

• 

• 

•r  P 

*3* 

CM 

CM 

03 

LT> 

CO 

oo 

o 

o 

PM 

03 

CM 

03 

i— 

in 

oo 

oo 

in 

CD  O 

*3* 

03 

03 

oo 

03 

CM 

00 

in 

00 

03 

03 

CO 

IM 

CO 

ro 

CM 

ro 

IM 

IM 

co 

4- 

03 

03 

03 

«3- 

«3- 

IM 

03 

CO 

in 

r— 

03 

IM 

*3* 

Mag 

(mb 

• 

*3- 

*3- 

• 

*3* 

LO 

1 

*3* 

1 

^3- 

*3* 

• 

in 

• 

*1* 

^3- 

*3- 

«3* 

• 

1 

• 

• 

*3- 

-C 

05 

cc 

CC 

CC 

oo 

4-> 

O 

*3- 

ro 

F 

X) 

1 

CO 

i 

03 

in 

CM 

ro 

ro 

no 

ro 

03 

1 

r— 

in 

Q.  E 

co 

ro 

s- 

CM 

ro 

in 

CM 

ro 

ro 

bo 

ro 

r— 

r— 

ro 

QJ  ^ 

O 

r— 

O'— 

z 

a> 

• 

QJ 

+-> 

03  S- 

C 

c 

1/3 

r— 

• 

• 

4-> 

c 

<u  0) 

03 

03 

u 

r 

• 

►— « 

■r 

1/3 

1/3 

o 

4- 

05 

S-  -*-> 

£X 

CL 

•f— 

• 

l/I 

in 

-C 

M 

NH 

n- 

O 

a 

03  C 

03 

03 

03 

■a  « 

M 

N 

o 

QJ 

L- 

c 

05 

• 

QJ 

*"D 

-3 

c 

0> 

03 

c 

3 

r— 

c 

1/3 

fO 

o 

■M 

1/3 

•—  o 

C.  03 

03 

LO 

03 

c 

4- 

u 

s- 

c 

•r— 

QJ 

O 

x: 

•r 

c 

o 

1/3 

HH 

J-  -r- 

• 

» 

c 

3 

r— 

03 

c 

QJ 

O 

s- 

-C 

>» 

O 

o 

CD 

L- 

XJ 

05 

* 

CL 

5-i 

3 

3 

c 

03 

cn 

03 

tn 

1/3 

•r 

“O 

QJ 

<_) 

c 

lO 

QJ 

QJ 

03 

•r 

O 

3 

JC 

C 

QJ  <D 

-C 

.c 

QJ 

03 

+-> 

c 

JC 

>— « 

+-> 

-C 

F 

«T3 

-C 

X 

o 

a. 

C  Oi 

jr  ^ 

05 

O 

-C 

05 

o 

C 

-M  O 

CO 

CO 

1 

03 

1 0 

•r> 

+-> 

in 

3 

4-> 

P 

4-» 

4-> 

4-> 

QJ 

•p* 

as 

QJ 

•M 

1/3 

1/3 

•r 

•p 

QJ  4- 

c5  o 

l-  r— 

3 

3 

X)  X) 

C 

o 

c 

C 

03 

4-» 

QJ 

3 

c 

L- 

C 

L. 

1/3 

03 

C 

03 

r— 

OJ 

• 

3 

05 

4- 

C 

XJ 

CD 

o  o 

>1 

>» 

•r- 

o 

3 

QJ 

<u 

r— 

03 

O 

OJ 

QJ 

<T3 

O 

QJ 

TO 

QJ 

3 

3 

1/3 

O 

r“ 

O 

4- 

o 

O 

QJ 

zo 

X  OH 

>- 

a. 

o 

< 

Oi  C 

oo  a. 

OO 

z 

3  Qc 

-d  a: 

cr>-> 

CO  < 

oc 

o 

z 

y 

Q1 

0) 

■o 

3 

LU 

LU 

3 

u 

3 

LU 

3 

3 

LU 

LU 

3 

O 

CO 

LO 

S 

LO 

03 

o 

CM 

3 

3 

3 

3 

3 

3 

<3- 

LU 

3 

4-» 

03 

CO 

CO 

»— 

03 

03 

03 

CM 

co 

CO 

«3* 

CO 

m 

oo 

03 

LO 

oo 

CM 

ct> 

CM 

O 

o 

CM 

o 

00 

03 

CM 

in 

03 

in 

in 

r— 

03 

CO 

CM 

ro 

CM 

c 

im 

ro 

ro 

CM 

CM 

oo 

*3- 

r*. 

oo 

co 

CO 

oo 

03 

CM 

<3* 

«3- 

in 

o 

» — 

F“ 

r“ 

r* 

—1 

0) 

“O 

z 

z 

z 

00 

z 

z 

z 

z 

z 

00 

oo 

Z 

3 

CM 

in 

3 

CM 

r>H 

o 

03 

^3- 

3 

OO 

LO 

LO 

z 

Z 

PM 

z 

z 

+-> 

oo 

*— 

•— 

03 

03 

o 

*— 

CO 

CO 

CM 

LO 

CM 

LO 

LO 

o 

n 

ro 

CO 

4-» 

CO 

CM 

CM 

d 

M 

00 

*3- 

ro 

CO 

_ _ 

o 

p_ 

r— 

r- 

o 

_ 

— 

CM 

o 

CO 

03 

ro 

ro 

ro 

CO 

in 

ro 

r— 

CM 

ct 

IM 

m 

CO 

*3- 

in 

—1 

8 

oo 

O 

im 

CM 

CO 

o 

00 

in 

o 

m 

o 

o 

r-  *r“ 

• 

• 

■ 

• 

• 

• 

• 

• 

• 

• 

• 

• 

■ 

• 

■ 

• 

03  4-> 

03 

03 

*3* 

m 

03 

00 

03 

o 

in 

CM 

03 

CO 

CO 

CO 

CO 

C 

u>  c 

o 

*3- 

O 

7; 

n 

CM 

*3" 

ro 

ro 

*3- 

7 

CM 

in 

T. 

t: 

CM 

CM 

o 

in 

ro 

in 

•r  *r 

03 

03 

p— 

03 

03 

CO 

in 

o 

00 

oo 

in 

00 

lO 

F— 

ao 

IM 

co 

r“ 

<✓>  CT> 

LO 

CM 

ro 

•— 

ro 

*7; 

•“ 

in 

ro 

in 

LO 

LO 

ro 

ro 

in 

o 

&- 

CO 

ID 

CO 

^3- 

M 

CM 

O'* 

ro 

CO 

o 

r- 

r- 

ro 

<3* 

00 

ro 

— 

co 

in 

O 

O 

o 

o 

O 

r— 

o 

O 

o 

&- 

&- 

i- 

s- 

i- 

&- 

S- 

s- 

i- 

i- 

C 

C 

c 

c 

c 

c 

c 

c 

c 

QJ 

£ 

£ 

03 

2 

03 

CL 

o. 

a. 

CL 

a 

3 

3 

3 

3 

3 

3 

3 

3 

3 

4J 

X 

X 

X 

< 

< 

«c 

< 

< 

rj 

ro 

X) 

>“D 

XJ 

rO 

"O 

•"D 

•X) 

03 

a 

ro 

CO 

DO 

o 

in 

in 

r— 

r^. 

r— 

— 

ro 

ro 

no 

n 

CO 

CO 

ro 

CM 

CM 

CM 

CM 

CM 

ro 

o 

o 

CM 

CM 

r* 

CM 

QJ 

03  r— 

* 

♦ 

« 

+->  Q. 

r— 

CM 

CO 

^3- 

n 

CO 

r^» 

00 

03 

O 

r— 

CM 

no 

*3- 

in 

ID 

PM 

GO 

03 

o 

<3  i 

r— 

” 

’ 

CM 

«/) 


-12- 


s 


TR  68-47 


V 


♦Three  station  location  using  arrival  times  at  TFSO,  UBSO,  WMSO. 


RELATIVE  AMPLITUDE 


Figure  9.  Filter  amplitude  responses 
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Figure  10.  Calcomp  plot  of  the  filtered  and  aligned  traces  for  Z1  through  Z10  from  data  sample 


An  arbitary  time-window  (generally  no  more  than  8  seconds)  was  established 
and  a  reference  time  selected  so  that  the  time  window  spanned  the  start  of  the 
signal.  In  this  window  the  largest  (peak-to-peak)/2  signal  amplitude  was 
determined  for  each  individual  element  and  the  phased  sum.  In  general,  this 
measurement  was  made  on  the  same  half- cycle  of  the  signal  on  each  trace. 

The  noise  analysis  was  accomplished  by  computing  the  rms  and  the  center-to- 
peak  range  of  the  noise  for  the  100  seconds  preceding  the  start  of  the  signal. 

Figure  11  shows  the  major  portion  of  the  TFOSAN  printer  output  for  PS3  of  data 
sample  6.  The  trace  amplitude  values  in  columns  2,  3,  and  4  are  given  in 
millimicrons.  Columns  5  and  6  are  signal-to-noise  ratios  and  column  8  contains 
the  delay  times,  in  seconds,  used  to  form  the  phased  sum.  For  each  statistic 
computed  in  the  analysis  of  the  individual  elements,  the  mean,  the  variance, 
and  the  standard  deviation  of  that  statistic  is  computed.  Program  TFOSAN 
analyzes  the  phased  sum  in  the  same  manner  as  ai  individual  element  is  analyzed 
and  compares  the  analysis  to  the  average  analysis  of  the  individual  elements. 
The  comparison  of  PS3  to  the  average  individual  element  for  data  sample  6 
shows  0.19  dB  signal  loss,  13.02  dB  rms  noise  reduction,  and  12.81  dB  signal- 
to-rms  noise  ratio  improvement. 

A  high-magnification  plot  of  each  phased  sum  was  obtained.  Figure  12  shows  a 
portion  of  PS3  for  data  sample  6.  The  magnification  of  the  plot  of  PS3  is 
approximately  eight  times  that  of  the  individual  traces  in  figure  10.  It  is 
interesting  to  note  that  PS3  snows  clearly  the  first  cycle  of  the  signal,  which 
is  completely  obscured  on  the  unprocessed  seismograms  and  only  faintly  visible 
on  the  filtered  individual  seismograms.  Figure  12  also  shows  the  high- 
magnification  plots  of  PS1  and  PS2  for  data  sample  6.  The  noise  reduction 
provided  by  PS2  over  PS1  and  PS3  over  PS2  can  be  seen  clearly  from  these  data. 

The  spectra  of  the  signal  and  the  noise  were  computed  for  each  of  the  elements 
and  averaged  over  the  data  sample.  A  maximum  lag  of  3  seconds  was  used  in  cal¬ 
culating  both  the  signal  spectra  and  the  noise  spectra,  resulting  in  estimates 
of  signal  power  and  no^se  power  at  the  same  frequencies.  The  spectra  for  the 
three  phased  sums  were  computed  in  the  same  manner,  except  for  the  averaging. 
The  spectra  for  the  average  individual  element  (IND) ,  the  three  phased  sums, 
and  the  filtered  summation  of  the  cross-linear  array  (CLSF)  are  shown  in 
appendix  5.  The  secondary  peak  in  the  noise  spectra  at  about  2.0  cps  was  a 
common  feature  of  all  of  the  noise  spectra.  Another  common  characteristic  of 
the  noise  was  that  the  average  peak  power  for  the  individual  elements  was 
generally  about  0.1  my^/cps.  For  the  20  signals  in  the  data  ensemble,  the 
signal  power  peaked  below  1.0  cps  on  18  signals,  with  the  majority  peaking  at 
0.67  cps.  It  should  be  noted  that  the  frequency  resolution  of  the  spectra  was 
about  0,16  cps,  which  allows  only  a  rough  estimation  of  peak  frequencies. 

The  complete  processing  scheme  for  a  data  sample  is  summarized  in  figure  13. 
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4.  DATA  ANALYSIS 


4.1  NOISE  REDUCTION 

For  each  data  sample,  the  noise  power  spectra  were  used  to  form  the  following 
ratios:  PS1/IND,  PS2/IND,  PS3/IND,  and  CLSF/IND.  These  ratios  indicate,  as 
a  function  of  frequency,  the  amount  of  noise  reduction  achieved  by  each  system. 
The  noise  reduction  ratios  for  each  data  sample  are  shown  in  appendix  1.  It 
is  evident  that  the  three  phased  sums  provide  essentially  the  vfT  noise 
reduction  at  the  higher  frequencies.  The  result  of  the  combination  of  the 
analog  and  the  digital  filtering  is  evident  in  the  noise  ratio  CLSF/IND, 
especially  at  the  higher  frequencies. 

The  noise  ratios  for  a  given  system  were  averaged  over  the  data  ensemble  and 
figure  14  shows  the  results.  From  the  average  noise  reduction  ratios,  we 
conclude : 


a.  Below  about  1.2  cps,  PS1  provides  about  3  dB  more  noise  reduction 
that  does  the  unfiltered  summation  of  the  cross-linear  array  (CLS) .  Above 
1.2  cps,  the  two  systems  are  about  equal. 

b.  PS2  provides  about  4  dB  more  noise  reduction  than  PS1  at  every 
frequency. 

c.  PS3  provides  about  2  dB  more  noise  reduction  than  PS2  at  every 
frequency. 

d.  The  noise  reduction  provided  by  the  CLSF  varies  considerably  due  to 
the  analog  filtering.  At  the  lowest  frequency  considered,  0.33  cps,  the  noise 
reduction  achieved  b>  the  CLSF  is  equivalent  to  that  of  PS2,  but  the  effective¬ 
ness  of  the  CLSF  decreases  with  increasing  frequency  until  at  1.0  cps  it  is 
equivalent  to  the  CLS.  For  frequencies  above  1.0  cps,  the  CLSF  provides 
increasing  noise  reduction  with  increasing  frequency.  At  the  highest  frequency 
considered,  3.0  cps,  the  CLSF  is  about  4  dB  more  effective  than  PS3.  Table  2 
summarizes  the  average  noise  reduction  provided  by  each  system  relative  to 

the  average  individual  element. 


Table  2.  Average  noise  reduction  at  3  frequencies 


•vsFreq 

Sys'^^ 

0.5  cps 

1.0  cps 

2.0  cps 

PS1 

5.1  dB 

7.0  dB 

9.7  dB 

PS2 

8.7  dB 

11.1  dB 

13.3  dB 

PS3 

11.3  dB 

13.7  dB 

15.0  dB 

CLS 

2.0  dB 

3.7  dB 

10.6  dB 

CLSF 

6.8  dB 

4.6  dB 

14.9  dB 
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The  rms  measurements  of  the  noise  on  the  i  , 

CLS,  and  CLSF  were  used  to  obtain  an  rms  noise  redlctivT^atio^*  PhaSwd  SUI"S’ 
relative  to  the  average  individual  element  Tahif f  V  f°r  each  system 
rms  noise  reduction  for  each  system?  glves  the  Value  of  the 


Table  3.  Rms  noise  reduction  ratios 


System 


PS1 

PS2 

PS3 

CLS 

CLSF 


Rms  noise 
reduction 

5.64  dB 
9.15  dB 
11.84  dB 
2.04  dB 
6.36  dB 


power°at'the  TolTr  'rLT'lt  ™  by  the  noise 

agreement  with  the  noiL^cS  X™d“?£ 

4.2  SIGNAL  LOSS 

signal^ower^pectr^were'use^to^et"  °f  the."°iSe  redUCti°"  rati°a-  «<• 

of  frequency.  Figure  15  shows  the  e™111®  s*8nal  l0SS  ratios  as  a  function 
These  curves  show  that  for  data  sLnle  6  th  werei  comPuted  for  data  sample  6. 

peak  power  was  1.4  dB  for  PS1  2  0  dB  for  pq?  ^o^n1?55  at  the  frecluency  of 
CLSF.  From  fWe  10  it  i«  LJa  ?  5  PS2'  1,9  dB  for  PS3>  «» d  3.6  dB  for 

on  the  individual  elements  are  the  ^  Sllght  differences  in  the  signal  coda 
phased  sums.  The  larger Lnal  r5“0n  f°r  the  si*nal  loss  on  the 

analog  filter.  ?he  sCnaHoss  °V  dUe  t0  the  r«P°nse  of  the 

appendix  2.  8  ratlos  for  each  data  samPle  are  shown  in 

(see  figure1?^). raThraverage8si?naiyitem  ""T®  avera8ed  over  the  data  ensemble 

range *of  iXT?rl^ 

Cbs  and  asAs  tLt^f  £ 

indicated^by'the^sipnal- loss 'ratios  “p*  Si8nalTJ0SS  ia  «  great  as  that 

over  the  data  s^l^L”  Te  XnTL £7“ 
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PS1 

0.37 

dB 

PS2 

0.55 

dB 

PS3 

0.62 

dB 

The  CLS  and  CLSF  were  not  included  in  table  4  because  the  amplitude  measurements 
of  the  signal  for  these  systems  were  calculated  using  Program  SPECSUB2  which 
calculates  the  maximum  center-to-peak  amplitude  in  the  signal  coda;  whereas, 
TFOSAN  calculates  the  maximum  peak-to-peak  amplitude  within  a  restricted  window. 
Because  of  this,  comparison  of  CLS  and  CLSF  signal  amplitudes  with  the  average 
amplitude  recorded  by  the  individual  elements  was  not  necessarily  meaningful  in 
all  instances. 

'jt 

4.3  SIGNAL-TO-NOISE  RATIOS 

The  signal  and  the  noise  power  spectra  for  each  system  were  used  to  obtain 
S/N  as  a  function  of  frequency.  The  curves  for  data  sample  6  are  shown  in 
figure  17.  It  should  be  noted  that  the  maximum  S/N  does  not  occur  at  the 
same  frequency  as  the  maximum  signal  power.  In  general,  the  S/N  curves  peaked 
at  about  1.17  or  1.33  cps.  The  S/N  curves  for  a  given  system  were  not  averaged 
over  the  ensemble  because  they  do  not  have  a  common  reference  point,  and  thus 
vary  from  data  sample  to  data  sample.  The  signal-to-noise  ratios  for  each 
data  sample  are  shown  in  appendix  3. 

Program  TFOSAN  computed  a  signal-to-noise  ratio  for  the  individual  elements 
and  the  phased  sums  by  taking  the  ratio  of  the  (peak-to-peak)/2  signal  amplitude 
to  the  rms  noise  value.  Table  5  gives  the  S/N  for  data  sample  6. 


Table  5.  S/N  from  Program  TFOSAN  for  data  sample  6 


System 


IND 

16.4 

dB 

PS1 

21.2 

dB 

PS2 

26.2 

dB 

PS3 

29.2 

dB 

4.4  S/N  IMPROVEMENT 

S/N  improvement  for  the  phased  sums  and  CLSF  was  computed  as  a  function  of 
frequency  by  taking  the  ratio  of  the  S/N  of  the  system  under  consideration 
to  the  S/N  of  an  average  individual  element  (IND) .  Figure  6  in  appendix  4 
shows  the  results  for  data  sample  6.  At  0.67  cps,  the  frequency  of  maximum 
signal  power,  the  S/N  improvement  is  4.0  dB  for  PS1,  7.8  dB  for  PS2,  10.8  dB 
for  PS3,  and  2.0  dB  for  CLSF.  At  the  frequency  where  the  maximum  S/N  occurs, 
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the  improvement  is  5.9  dB  for  PS1,  9.1  dB  for  PS2 
for  CLSF.  * 


11.5  dB  for  PS3,  and  2.0  dB 


CTsemble1"arHVfhent  ra^0S  f°r  !ach  system  Were  averaSed  °ver  the  data  sample 
ensemble,  and  the  results  are  shown  in  figure  18.  The  ratios  show  that  the 

PS2  Lki  1V£?  following  order,  from  poorest  to  best:  CLS,  CLSF,  PS1, 
“d  PS3-  Table  6  summarizes  the  S/N  improvement  provided  by  each  system 
relative  to  the  average  individual.  7  system 


Table  6.  S/N  improvement  at  3  frequencies 


System 

0.67 

cps 

PS1 

4.1 

dB 

PS2 

7.5 

dB 

PS3 

10.0 

dB 

CLS 

1.1 

dB 

CLSF 

2.8 

dB 

I ■ 0  cps  1.17  cps 

5.6  dB  6.4  dB 

9.3  dB  10.0  dB 

11.9  dB  12.4  dB 

1.7  dB  2.7  dB 

2.3  dB  2.9  dB 


It  should  be  noted  that  the  drop  in  the  S/N 
(^2.0  cps)  is  due  to  the  general  absence  of 


improvement  at  the  higher  frequencies 
signal  power  at  these  frequencies. 


Average  S/N  improvement  was  also  computed  from  the  S/N  improvement  factors 
°“  P.  by  Program  TFOSAN  for  the  three  phased  sums.  Phased-Sum  1  shows  an 
average  S/N  improvement  factor  of  5.26  dB,  PS2  shows  8.58  dB  improvement,  and 
PS3  shows  an  improvement  of  11.25  dB.  These  values  compare  favorably  with  the 
S/N  improvement  ratios  at  1.0  cps  (table  6)  X 


4.5  CONCLUSIONS 

The  following  conclusions  can  be  made  from  the  data  analysis: 

1.  The  three  phased  sums  and  CLS  suppress  the  noise  by  a  factor  of  VTT 
for  frequencies  greater  than  about  1.5  cps.  At  1.0  cps,  the  average  noise 
I®  nni°”  relatlve  t0  an  average  individual  element  is  7.0  dB  for  PS1  11  1  dB 
d“  nfor  pS3,  3.7  dB  for  CLS,  and  4.6  dB  for  CLSF.  For  the  pass 
n  ??  4  t0  3,0  cps*  the  average  rms  noise  reduction  is  5.64  dB  for  PS1 

9.15  dB  for  PS2,  11.84  dB  for  PS3,  2.04  dB  for  CLS,  and  6.36  dB  for  CLSF. 

i  n  ^ j  amplitude  loss  in  the  30-second  P-wave  coda  is  between 

l.u  and  2.0  dB  in  the  frequency  range  of  greatest  signal  power  (0.5  to  1.0  cps) 

r?Qpthe^hree  fha!67  the  CLS>  .  Due  t0  the  resP°nse  of  the  analog  filter, 

F  suffers  about  3.0  dB  signal  loss  in  the  same  frequency  band.  Using 

(peak-to-peak)/2  signal  amplitude  values  measured  on  a  dominant  half-cycle  in 
a  time  window,  the  average  signal  loss  relative  to  an  average  individual  is 
0.37  dB  for  PS1 ,  0.55  dB  for  PS2,  and  0.62  dB  for  PS3. 
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Figure  18.  Ensemble  averaged  S/N  improvement  ratios 


3.  From  a  S/N  improvement  standpoint,  the  five  systems  that  were  studied 
rank  as  follows,  from  poorest  to  best:  CLS,  CLSF,  PS1,  PS2,  PS3.  The  maximum 
S/N  improvement  generally  occurred  at  about  1.17  or  1.33  cps  on  the  three  phased 
sums.  At  1.17  cps,  the  average  S/N  improvement  relative  to  an  average  indivi¬ 
dual  is  6.4  dB  for  PS1,  10.0  dB  for  PS2 ,  12.4  dB  for  PS3,  2.7  dB  for  CLS,  and 
2.9  dB  for  CLSF.  For  S/N  computed  with  (peak-to-peak)/2  signal  amplitudes  and 
rms  noise  values,  the  average  S/N  improvement  is  5.26  dB  for  PS1,  8.58  dB  for 
PS2,  and  11.25  dB  for  PS3. 
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NOISE  REDUCTION  RATIOS  FOR  EACH  DATA  SAMPLE 
RELATIVE  TO  AN  AVERAGE  INDIVIDUAL  ELEMENT 
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Figure  10.  Noise  reduction  ratios  for  data  sample  10 
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Figure  12.  Noise  reduction  ratios  for  data  sample  12 


Figure  14.  Noise  reduction  ratios  for  data  sample  14 


Figure  17.  Noise  reduction  ratios  for  data  sample  17 
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igure  19.  Noise  reduction  ratios  for  data  sample  19 
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SIGNAL  LOSS  RATIOS  FOR  EACH  DATA  SAMPLE 
RELATIVE  TO  AN  AVERAGE  INDIVIDUAL  ELEMENT 
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Figure  4.  Signal  loss  ratios  for  data  sample 
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Figure  5.  Signal  loss  ratios  for  data  sample 


Figure  7.  Signal  loss  ratios  for  data  sample 
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APPENDIX  3  to  TECHNICAL  REPORT  NO.  68-47 
SIGNAL -TO -NOISE  RATIOS  FOR  EACH  DATA  SAMPLE 
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Figure  17.  Signal-to-noise  ratios  for  data  sample  17 
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SIGNAL-TO-NOISE  IMPROVEMENT  FOR  EACH 
DATA  SAMPLE  RELATIVE  TO  THE  SIGNAL-TO-NOISE  RATIO 
FOR  THE  AVERAGE  OF  THE  INDIVIDUAL  ELEMENTS 
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Figure  3.  Signal-to-noise  improvement  for  data  sample 
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Figure  6.  Signal-to-noise  improvement  for  data  sample 
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Figure  7.  Signal-to-noise  improvement  for  data  sampl 
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Figure  15.  Signal-to-noise  improvement  for  data  sample  15 
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Figure  17.  Signal-to-noise  improvement  for  data  sample  17 


Figure  18.  Signai-to-noise  improvement  for  dat 


TR  68-47,  app  4 


Figure  19.  Signal-to-noise  improvement  for  data  sample  19 
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NOISE  AND  SIGNAL  POWER  SPECTRA  FOR  EACH  DATA  SAMPLE 
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Figure  1.  Noise  reduction  ratios  for  data  sample 
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Figure  6.  Noise  reduction  ratios  for  data  sample 
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Figure  9.  Noise  reduction  ratios  for  data  sample 
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Figure  11.  Noise  reduction  ratios  for  data  sample  11 
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Figure  14.  Noise  reduction  ratios  for  data  sample  14 
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Figure  15.  Noise  reduction  ratios  for  data  sample  15 
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Figure  17.  Noise  reduction  ratios  for  data  sample  17 
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Figure  19.  Noise  reduction  ratios  for  data  sample  19 
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Figure  20.  Noise  reduction  ratios  for  data  sample  20 


